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Abstract 

The writhe of a space curve fragment is considered for various 
boundary conditions. An expression for the writhe as a function of 
arclength for an arbitrary space curve is obtained. The formula is built 
on the base of closing the tangent indicatrix with a geodesic. The cor- 
responding closure of a curve in 3-space is explicitly constructed. The 
addition rule for writhe is formulated. A relationship connecting the 
writhe with the GauB integral over the open curve is presented. The 
single and double regular helical shapes are examined as examples. 
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1 Introduction 

The term "writhe" (or "writhing number") was first proposed by Fuller 
for a quantity Wr that arises as a difference between the linking num- 
ber Ck and the twist (or "twisting number") Tw of a closed ribbon in the 
Calugareanu- White-Fuller formula ||2|, ^, ^ 

Wr = Ck- Tw (1) 

The sense of this very simple-looking and very famous (isn't the second a 
consequence of the first?) relationship is that the right-hand side, though 
defined for a ribbon, depends only on its central curve. 

Since Wr is a characteristic of spatial complexity of a curve, it makes 
this quantity worth to consider while examining lengthy physical objects. 
In particular, the values of Wr have been computed in a number of works 
for various models of large-scale structure of DNAs (e.g., |5|, ^]) as well 
as for experimental data on these molecules [^. RNA tertiary structures 
and the protein folding are other neighbour areas for an application of the 
geometrical and topological tools developed in the DNA studies How- 
ever, the application of the writhe is confined formally by smoothly closed 
shapes though quite a lot of interesting objects have their ends not joined 
(or joined non-smoothly) . One should mention here the studies based on 
the recently developed experimental techniques of manipulation with single 
DNA molecules |^. The modelling curves require an appropriately defined 
measure to characterize their arrangement in space. Therefore, it seems to 
be helpful to generalize the notion of the writhe onto non-closed curves and 
their fragments. Indeed, such trials (explicit and implicit) have been made 
[10, 11, 1^, |13|, 14, 15, 16, 0. In some cases the definitions of the writhe 



suggested in different works are not consistent to each other. 

The aim of this paper is to give a consistent and natural generalization 
of the notion of writhe to an arbitrary fragment of a curve and to derive 
explicit formulas for its computation as a function of the arclength. The 
basic idea consists in a construction of the closure of the fragment under 
consideration in such a way that it would correspond to the closure of the 
tangent indicatrix by an arc of a great circle, as it has been proposed by 



Maggs [15, |l^. In some sense, we can think then about the writhe as being 



locally defined (cf. M) 
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2 Notation and preliminaries 

We start our consideration with a smooth non-self-intersecting curve A = 
r{s) : [0,L] — > of class C^, s being the arclength. We assume that 
the segment has a natural orientation in the direction of the arc coordinate 
increase. In what follows, we shall denote a reverse-oriented segment by 
the minus sign: —A. The concatenation of two segments A and B having, 
respectively, the common ending and starting point will be written as A + B. 

A vector function u(s) : [0, L] ^ = {z G M'^; |z| = 1}, u € C, may be 
chosen such that s, u(s) • r'(s) = 0, Vs; ' denotes the derivative with respect 
to s. Let e > be small enough so that the ribbon = {r+/iu, —e<fi<e} 
does not cross itself. 



2.1 The twist, the writhe and the hnking number 

The twisting number (or the twist) of the ribbon (i.e., of the pair (r, u)) 
is defined by 



Tw(r,u) = — [ r' X udu. 







If r G and r" ^ 0, then the twist of a ribbon may be decomposed into 
the twist of the Frenet frame plus the twist of the ribbon relative to the 
Frenet frame [^]: 

L 

Tw{r, u) = TwF{r) + ^ ^ #• (2) 



The angle (f) = (j){s) is an angle between u and the principal normal. By the 
Frenet ribbon is meant a special one defined by the principal normal vector 
to the curve. The Frenet ribbon is defined uniquely by the space curve r(s), 
if r" 7^ 0. The twist of the Frenet ribbon is 

L 

1 



TwFir) = ^ I r{s)ds, (3) 







where r(s) is the torsion of the curve r(s). Clearly, for planar curves, 
TwF = 0. 
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Now consider the smoothly closed curve: r(0) = r(L),r'(0) = r'(L). The 
quantity called the writhe may be expressed as the double integral 



^ f ? {r{si)-r{s2))-{t{si ) X t(s2)) , 

4^7 7 |r(si)-r(s2) 




where t = r'(s) is the tangent vector and si,S2 are the arclengths. The 
right-hand side of Eq. (^) is the Gaufi linking integral in the singular case 
as being over all distinct pairs of points on one curve. The writhe depends 
exclusively on the shape of the curve. 

For two closed curves A and Ar\ B = 0, the Gaufi linking integral 
gives an integer-valued topological invariant 

Ck{A,B) = — / ■ — — — dsids2, (5) 

^■^ J J \^A{si) - rB[S2)r 

B A 

called the linking number. 
2.2 The basic relations 



Two theorems by Fuller [gy, 21 1 are valid. They provide the means to 
compute Wr efficiently. 

1. Let B = r{s) G be a closed oriented space curve with its tangent 
r^(s), s the arclength. The tangent traces out a closed curve B{s) on the 
unit sphere which is piecewise of class C"^. The curve B{s) is divided into a 
finite family of non-self-intersecting closed piecewise space curves. Each 
curve of this family then encloses a domain defined so that the geodesic 
normal points into its interior. Let Sb be the sum of the areas of these 
domains (the components are counted with multiplicity determined by how 
many times the corresponding domains are encircled by the curve). Then 

Sr 

Wr{r) = — - 1 mod 2. (6) 
2tt 

2. Let ro(^) and ri(^) be two closed non-self-intersecting curves in 

of class C^, 9 £ [0, 0] being the common regular parameter. Let there exist 
a continuous deformation F(A) : {9, A) T^xi9), A G [0, 1], such that 

1) F(0) = ro(e) andF(l) = ri(0), 

2) rx{9) G is a non-self-intersecting curve for all A G [0, 1], 
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3) the tangent tx{9) is continuous in A, 

4) |z(to(^), tA(^))| < vr, y{e, A) G [0, e] x [0, 1]. 

Then the difference of the writhes of the two curves is 



We should stress that both theorems are apphcable to closed curves. 
Nevertheless, the works occur (e.g., |2^ ]) in which the writhe of a linear 
(non-closed) curve is computed by means of Eq. (^) that leads generally to 
improper results. 

Let the ribbon (r, u) be also closed: u(0) = u(L). Denote by Ck{r — 
//u, r + //u) = £A:(r, u) the linking number of the two boundary curves 
r — jjLU and r + //u. For e small enough, Ck(r, u) does not depend on e. This 
justifies omitting e in the following. In other words, we shall be dealing with 
arbitrarily narrow ribbons. 

The famous Calugareanu- White-Fuller theorem [^, ^, |^ claims that the 
difference of the linking and twisting numbers is the writhe: 

Ck{r, u) - Tw{r, u) = >Vr(r). 



2.3 The Gaufi integral over self-intersecting and non-smooth 
curves 

We are going to show that the integral in Eq. (^ ) is well-defined also for 
curves that are smooth only piecewise. To demonstrate this, it is sufficient to 
consider a fragment of a curve in the vicinity of the point where the tangent 
vector is not continuous. If the double integral does not diverge in this point 
then, as a consequence, it also exists for curves crossing themselves in a finite 
number of points. 

As a general case, we take the model of two curves (which may be thought 
of as two smooth fragments of the same curve) ri(si) and r2(s2) that have 
one common point O : ri(0) = r2(0); si and S2 are the arclengths, Sj G 
[0,/j],i = 1,2. Generally, ti(0) ^ t2(0). Each of the two curves is of 
class and, near O, it may be locally described in its own Frenet frame 
(tio,njo,bio),i = 1,2, as 

r.(«.) = [s, - ^sl ^sl + ^sl ^.?) + 0(4), ^ = 1, 2, (8) 
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where Kj and Tj are the curvature and torsion of the i-th curve in the point 
O. (If the point O is an (isolated) inflection point, it is still possible to 
define the local frame as the limiting Frenet for Si 0+.) Differentiation 
of Eq. (H) yields 



Uis,) = (l - ^sl ^,s, + ^sl ^si^ + 0{s\\ i = l,2. (9) 

Let A = {ajk} be an orthogonal matrix of orientation of the second 
Frenet frame with respect to the first one: 





t20 


n2o b2o 


tio 






nio 






bio 







Then the integrand of the Gaufi integral (Eq. (|5D) 

T ( ^ (r-l(gl) -r2(g2)) • (tl(si) X t2(52)) 

|ri(si) - r2(s2)r 

may be represented in the form 

, , , (ri(si) - Ar2(s2)) • (ti(si) X ^t2(s2)) 

|ri(si) - Ar2{s2)r 
_ At2 ■ (ri X ti) +ti • A{r2 x t2) 
Iri-ArsP 

With the help of Eqs. (^), (|9|), it is possible to obtain an approximation for 
small si, S2 

IWr{si,S2) = |^(a3iK;iSi + ai3K2s|) + 
+ ^[(a3l4 - a2lKlTi)s? + (ai3K2 - ai2K2T2)s2] + 
+ ^KiK2(a32Sl + a23S2)siS2 + X 
x[si - 2aiiSiS2 + ^2 ~ («2l'«lSl + ai2K2S2)siS2 + ©4]"^, (10) 

where the terms of the 4-th order and higher are denoted by 



1=0 
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Consider first the case when afi 7^ 1. It is convenient to introduce the 
new variables p, cj) {p > 0, (p £ [0, |]) such that 

si = — p(Vl - ail cos 4> + a/1 + ail sin</)), 
S2 = — p(a/1 + ail sin (f) - a/I - an cos (/>). 
Then, with the use of Eq. ([lO|), we have 

Iwrisi,S2)dsids2 = j j -^>Vr(/5, 0) y'T^-ofi p dp d(j) = 



l-ai^J J p3[l + yif^{(j))p + 0{p^)]2 
As p — > 0, the integrand of the last integral approaches 



^o(0) = ^[(aai^i +ai3K2)(l - an cos 20) + y 1 - a2^(a3iKi - ai3K;2) sin 20]. 

It vanishes only if 031 ki = and ai3K2 = 0. This happens when either 1) the 
both curves have O as an inflection point or 2) one curve has an inflection 
point and its tangent lies in the tangent plane of the other one or 3) the 
tangent planes of both curves coincide. 

Now let On = —1. The meaning of this is that the curves ri and r2 
form a curve having a continuous tangent. We introduce the angle x which 
parametrizes the relative orientation of the Frenet frames at O (or one can 
interpret x as a discontinuity angle between principal normals). Then we 
have for the entries of the orientation matrix: 012 = a2i = ai3 = 031 = 0, 
a22 = —033 = cosx, a23 = a32 = sinx- Eq. ([lO|) takes the form 

^KiK2Smx{si + S2)siS2 + 04 

IWr[Sl,S2) = 3 . 

[(^1 + 52)^ + 04]-^ 

Let us redefine the variables p and (p: si = psincj), S2 = pcoscj) {p > 0, 
G [0, f]), then 

IWrisi,S2)dSidS2 = j j -^Wr(p, 0)p dp d(j) = 

\ Ki K2 (sin (f> + cos (f>) sin (f) cos cpp'^ + O {p^ ) 

3 dp dcj), 

p3[(sin (f) + cos 0)2 + ©(p^)] 2 
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from where we can see that the integral exists. 

The case an = 1 corresponds to the cusp points and requires a special 
consideration which is beyond the scope of the present paper. Therefore we 
exclude cusps from the following analysis. 



3 A periodic curve 

On the one hand, the treatment of the periodic curve case is the most resem- 
bling to that of the closed configuration. On the other hand, the periodic 
curves are often met in various physical and, in particular, biomechanical 
models. For example, the linear helices and interwound structures are com- 
mon in the DNA modelling (e.g., see [p3| ). 



3.1 A formula for the writhe 

Let rA{s) be periodic in space: r^(s + La) = t^al + ryi('S), La the period, 
rAL = const. For the tangent tA{s) = r^(s), it implies tyi(s + La) = 
t^(s). Denote t^o = t^(0) and t^i = tyi(Lyi). Consider the segment 
^ : < s < La- The tangent indicatrix is a closed curve which sweeps 
out an area Sa on the unit sphere S'^ (in the same sense as in the Fuller 
first theorem). Now construct a ribbon Ra for the segment A in the same 



manner as it is done in the proof of the Fuller first theorem [ 21 1 . Generally, 
the ribbon Ra is determined by the unit principal normal vector n^(s). 
However, we should account for the possible inflection points of A with the 
discontinuous normals. Because of this, in the vicinities of the inflection 
points the ribbon Ra is to be arbitrary modified to make it continuous 
with a new modified generating vector u^. Due to the periodicity property, 
u^o = ua(0) = ua{La) = uai- For the sake of simplicity, we shall be also 
assuming, here and further, that u^o = n^(0) and u^i = n^(-L^). 

The twist of Ra, which is well defined for non-closed ribbons as well, 
satisfies the equation 

Tw{rA,UA) + ^= mod 1 (11) 
zvr 

(the above follows from the proof of the Fuller first theorem which, in turn, 
is a consequence of the Gaufi-Bonnet theorem). 

At each point of the ribbon A we can define an orthonormal frame 
{t^(s), u^(s), tyi(s) X Uyi(s)}, which is an element of SO{3). Thus, the 
ribbon A corresponds to a closed loop A in the space of orientations SO (3) 
[p4|. Following Hannay p^, we define an integer number G Z2 which is 
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zero if A may be continuously deformed to a point (i.e., if it is contractible) 
and = 1 otherwise. Then Eq. ( |TT| ) is a consequence of the Hannay 
formula 

'7ui(r^, u^) + — = mod 2. 

zvr 

Now build a planar non-self-intersecting segment C : rc(s) G C^, for 
< s < Lq, that closes the segment A. That means that rc{0) = ryi(Lyi), 
rc{Lc) = TAiO), tc(0) = tAi, tc{Lc) = tAo (tc = r'cis)) and C lies in the 
plane determined by t^o aiid r^L (if ^Ao II r^/,, then we choose the plane 
spanned by t^o and t'^g). The closing segment C may be chosen such that 
its e- vicinity does not cross the ribbon Ra- In addition, we require that the 
principal normal at the beginning and at the end points be well defined and 
have the same direction: nc(0) = nc(Lc). Then the number of inflection 
points of C is even. 

There exists a ribbon (rc(s), vc(s)) lying entirely in the plane of C 
with its ends defined by the normal vectors vc(0) = nc(0) and vc(Lc) = 
nc(Lc). Clearly, its twist is zero. (In the case of no inflection points, this 
ribbon coincides with the Frenet ribbon.) Construct another ribbon on the 
curve C with a vector uc(s), uc(s), uc(s) • tc(s) = 0, which differs from 
the vector vc(s) by a constant angle: uc{s) • vc(s) = cos(/)o = const, and 
uc(0) = u^i- It is evident that uc{Lc) = uc{0), uc{Lc) = u^o and 
Tw{rc,uc) = 0. 

Thus, joining the (oriented) ribbons A : (r^jU^) and C : {rc,uc) gives 
a continuous closed ribbon which we shall denote hy A + C. The linking 
number of a closed ribbon is an integer and, by applying the Calugareanu- 
White-Fuller relation to ^ + C, we have 

CkA+c = WrA+c + '^w{rA,UA) +Tw{rc,uc). (12) 

Prom Eq. (12) we obtain Wva+c = ^kA+c — '^w{rA,'ViA) and, taking 



WrA+c = ^ mod 1. (13) 
zvr 



into account Eq. (O 



In 1 26] it was shown that the last relation may be refined as follows 



WrA+c = ^ + Tnc + l mod 2, (14) 

ZTT 



where Trie is the turning number of the closing curve. 



E.L.Starostin. On the writhe of non-closed curves 



11 



Since The is integer, the fractional part of the writhe does not depend 
on the closing curve C and we may write 

WrA = 7— mod 1. 
zvr 

We see that in the case of a periodic curve the same formula is valid for 
the writhe of a period as in the case of a smoothly closed loop. 

Note that the tangent indicatrix of the joined ribbon differs from that of 
A by an appendage corresponding to the planar segment C. This appendage 
has an area of 27r The, hence, it does not contribute to the total spherical 
area of the domain enclosed by A + C, counted modulo 2tt. In other words, 
Sa+c = Sa mod 27r. 

If the writhe of a period is to be calculated by using the double integral 
formula of Eq. (Q), then at first the segment of the period should be closed by 
an additional planar segment as described above and after that the formula 
Eq. (^) should be applied to the entire closed loop, though the double 
integral in Eq. is also well defined for open segments. 

Another way to compute the value of the writhe WrA is the application 
of Eq. (0): 

WrA = CkA+c - TwivA, ua), (15) 

which is valid for any appropriate closing segment built as described above. 
Note also that Eq. (^) is an exact formula. With the help of the relation 
between the linking number and the contractibility number |^6|, Eq. (|l5| ) 
may be written in the form 

La 

Wta = N + Tnc + 1 / r'^ X UA • duA mod 2. (16) 

27r J 



Eq. (16) may be useful since one often needs to know only the fractional 



part of the writhe. 
3.2 Example 

Consider a 27r-periodic spatial curve 



1 si 
r = (— cos s, — - cos 2s, - + - sin 2s) (17) 



(Fig. |l|). Its tangent vector is 



t = (sin s, - sin2s, cos^ s) (18) 
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and the tangent indicatrix is a 8-shaped self-intersecting curve consisting of 
the two identical loops (Fig. ^). Let 5i be the spherical area inside one loop, 
then the complementary area of the other loop equals 47r — 5i. Thus, the 
sum of the areas of the two domains defined as in the Fuller first theorem, 
amounts to An, and, according to Eq. (|l4[), Wr = mod 2. 




Figure 1: One period of the curve of Eq. (p!7|) (solid line), its planar closure 
(dashed line) and the second edge of the Frenet ribbon (thick line). 

The last conclusion is evidently true for every periodic curve, the tangent 
indicatrix of which takes a symmetric form of the figure-8 with the identical 
loops. (In fact, the loops need only have equal areas.) 

This result may be easily verified (and refined) by applying Eq. ([TsD. To 
this end, we first construct a Frenet ribbon with the principal normal vector 

1 dt 1 , 
n = — — = —(cos s, cos 2s, — sin 2s), 
Kas K 

the curvature k = Vl + cos^ s. One edge of the Frenet ribbon r + en is 
shown in Fig. || as a thick line. The one-period segment may be closed with 
a planar curve as described above and an additional segment of the ribbon 
is to be built to close the Frenet ribbon of the initial curve. It is easy to 
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Figure 2: The tangent indicatrix (Eq. (|T8|)). 

check that the linking of the two edges for the entire closed loop is zero. Now 
compute the twist of the open Frenet ribbon by using Eq. (|3|) into which the 
torsion r = ■\r'r"r"' = — sins '^"'^^ ^"j"? is substituted. Clearly, 

COS^ S + 1 ' 

TwFir) = Tis)ds = 0, 



and by Eq. (jl^) we obtain finally Wr = 0. 

4 An open segment of a curve with the parallel 
initial and terminal tangents 

We again deal with the segment A of a spatial curve rA(s),0 < s < La- 
We now make no assumptions on the behaviour of this curve outside the 
segment A. Only one additional condition is to be satisfied: the tangent 
vector at the beginning is the same as one at the ending point: t^o = 
t^(0) = tA^LA) = tAi- Then the tangent indicatrix is again a closed curve 
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sweeping out the area Sa- We construct the ribbon Ra for the segment 
A by using the unit vector u^(s), Uyi(s) • t^(s) = 0, in the same manner 
as in the case of the periodic curve. The difference now is that, generally, 
u^o = u^(0) 7^ ua{La) = UAi- We denote by 7 the angle from u^i to u^o 
(actually, from uai to n^o) (Fig- §)• Note that both vectors lie in the plane 
orthogonal to t^o = tAi- 




Figure 3: The closed tangent indicatrix A of the segment A. 



Due to the GauB-Bonnet formula |27|, 



Sa + 27rTw{rA,UA) +7 = mod 2tt. (19) 

The next step is to build a planar non-self-intersecting segment C : 
rc{s) G C^, for < s < Lc, to close the segment A. C lies in the plane 
spanned by t^o and val (or t'^(O), if t^o || yal) and rc(0) = rA{LA), 
rc{Lc) = r^(0), tc(0) = t^i, tc(Lc) = t^o- We can also require that 
C have an even number of inflection points, all of them being interior. In 
addition, C is such that it is separated from the ribbon Ra by the distance 
more than e. 

A ribbon on C exists that lies in the same plane as C and with its 
end vectors coinciding with nc(0) = nc{Lc)- Denote its generating vector 



E.L.Starostin. On the writhe of non-closed curves 



15 



by vc(s). We build a ribbon on the curve C by using a vector uc(s), 
uc(s) ■'tc('S) = 0, which is the vector vc(s) turned through a constant angle 
(pQ-. nc{s) ■ vc(s) = cos(/)o = const. Choose (po such that uc(0) = uai- By 
the construction of C, uc{Lc) = uc(0) and the ribbon C has zero twist: 
Tw{rc,uc) = 0. 

We see that the joined ribbon A + C may have only one point of discon- 
tinuity, namely, the starting point of A (corresponding to the terminal point 
of C). At this point the twist angle between \ic{Lc) and u^o is exactly 7. 

The following relation is valid, according to [^8|, for the discontinuous 
ribbon A + C (see also Appendix A for a different proof) 

2tt{Twa+c + WrA+c) +7 = mod Itt. (20) 

The sum Twa+c + VVr^+C) which is not integer, may be interpreted as 
the linking number of a cord. The latter is a generalization of a notion of 
the ribbon and was introduced by Fuller in |2C]. 



Since the twist is additive, Twa+c = Twa + Twc, Twa = T'w{rA,UA), 
Twc = Tw(rc,uc) = 0. 

Combining Eq. (Il^) and Eq. (|ol) leads to 



WrA+c = ^ mod 1. (21) 

By the argument similar to that used in the case of the periodic curve it 
may be shown that 

WrA+c = ^+Tnc + l mod 2. (22) 
zvr 

Again, the tangent indicatrix of C does not change the area of the domain 
enclosed by A + C, counted mod 27r, and we conclude that the fractional 
part of the writhe does not depend on the particular shape of the closing 
curve constructed under the above conditions: 

WrA = ^ mod 1. (23) 

The fractional component of Wr may be also computed by means of 
Eq. (IgD: 

WrA = -TwA - ^ mod 1. (24) 
27r 

Note that Eq. ([2^ ) does not refer to the closing segment. All what is 
necessary to know is the twist for the original segment A and the angle 7 
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between the normals n^o SLud n^i. An analogue to Eq. was used in 
the analysis of the elongation of a supercoiled DNA molecule carried out 
by Bouchiat and Mezard |jl^, Q (though their angle x is measured in the 
opposite direction to 7). 



5 An arbitrary open segment 
5.1 The general case 

The most general case takes place if we dismiss the equality condition of the 
tangent vectors at both ends of the spatial curve segment A : rA{s),0 < s < 
La- In other words, the tangent indicatrix has not to be closed any more. 
Nevertheless, we can construct the ribbon Ra '■ (r^, u^) in exactly the same 
manner as we did it previously. The orientations of the generating vector 
at the ends may be different. 

In order to get a measure for the writhe of A consistent with the above 
considered cases we choose to close the tangent indicatrix with a geodesic 
[p^ , |l^ ] . This choice is natural and it is supported by treatment of analogous 
problems in optics and quantum mechanics [^]. (In the generic case t^o 7^ 
ibt^i, there are two possible geodesies, we take one of them; the case t^o = 
—tAi will be discussed later.) 

Let G be a planar segment rG(s),0 < s < Lq-, in the plane determined 
by tAo,tAi such that rG(0) = va{La) and r'g(O) = t^i, '^'(.{Lg) = t^o- We 
may always require that G have no inflection points. 

We build a ribbon Rq based on the segment G and the vector ug'(s), < 
s < Lg, which is the principal normal nG{s) turned around the tangent 
through the constant angle —71 such that ugo = ug'(O) = u^i. Therefore, 
the two ribbons Ra and Rg are continuously glued to produce the joined 
ribbon Ra+g- 

We now can see that the assumption of the previous case is satisfied for 
the joined curve A + G: it has the same tangents at the ends. By applying 



Eq. (ID to A + G, we obtain 



WrA+G = ^^+'rnc + l mod 2. (25) 
zvr 



It should be noted that although Eq. (|25| ) is valid mod 2, actually 
we may define and compute only fractional part of the writhe because of an 
arbitrariness of one of two closing geodesies chosen. To put this another way, 
we are able to determine the area only mod 27r. Therefore, the writhe of 
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an open segment may be determined by the following relation: 

>Vr^ = ^d±£ mod 1. (26) 



27r 



Instead of Eq. (24), we come to 

T 

WrA+G = -TWA+G - TT 1' (27) 

zvr 

where Twa+g = Twa + Twg- But Twg = Tw{yg-,'^g) = by the same 
argument as for the closing segment for a periodic curve. 

The angle 7 is the angle measured from mgi = vlg{Lg) to u^o- It may 
be represented as the sum 7 = 71+70, where 70 is an angle from the normal 
nci to u^o(= n^o) (Fig. ||). 




Figure 4: The tangent indicatrix A of the segment A closed with a 
geodesic G. 

Finally, we arrive at 

WrA = -TwA - ^ ^° mod 1, (28) 
27r 

where the angles 71 and 70 are determined by the equations 
cos 71 = n^i • nco, sin 71 = (hai x nco) • tAi, 
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cos 70 = n^o • nci, sin 70 = (nci x hao) • tAo- 

The normal vectors nco and nci may be easily expressed as functions 
of the initial and terminal tangents of A. Indeed, the vector nco lies in 
the plane spanned by the vectors tAo and tAi- Also, nco • tAi = and 
nco • tAo > 0. Hence, 

tAO - (tAO ■ ■tAl)tAl 
- (tAO • tAl)2 

Similarly, the vector nci lies in the same plane and nci ■ tAo = and 
nci • tAi < 0. Therefore, 

tAl - (tAO • tAl)tAO 
Vl-(tA0-tAl)2 

The angles 71 and 70 can be found from their trigonometrical functions 

nAl • tAO . bAl • tAO 

COS71 = — , sm7i = — , (29) 

v/l-(tA0-tAl)2 Vl-(tA0-tAl)2 

-HAO • tAl . bAo • tAl 

COS70 = — , sin7o = — , (30) 

Vl-(tA0-tAl)2 Vl-(tA0-tAl)2 

where bA(s) = tA(s) x nA(s) is the binormal vector and bAo = bA(0), 
bAl = bA(LA). 

We can conclude that Eq. (|28|), together with Eqs. (p9[), (pO[), provides a 
means to compute the fractional part of the writhe for an (almost) arbitrary 
curve with open ends. 

Remark 1 . It follows from the above that the writhe of a curve segment 
(closed or open) whose tangent indicatrix is geodesic and such that tAo + 
tAl 7^ 0, is an integer. In other words, the writhe of any planar curve is 
always integer. 

Consider now the case tAo = — tAi- The closing geodesic G is not deter- 
mined uniquely then. If we examine the behaviour of the tangent indicatrix 
in the vicinities s = and s = La (also paying attention to the neighbour 
curves on close to the critical one), then we see that the plane of the clos- 
ing geodesies may rotate through ~ vr as the length of the segment changes 
so that the critical point tAo = — tAi is passed. In mechanics terms, this 
phenomenon may be called flipping of the closing segment. It is the critical 
point where the choice of one geodesies based on the continuity argument is 
no more valid. 
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Speaking more strictly, the value of the writhe for a segment with the 
oppositely directed ends is not determined. It may be ascribed a value which 
is average of the two limits taken as the length of the segment is pre-critical 
and post-critical. That means that the great circle plane for the closing 
geodesies in the critical point has to be taken orthogonal to the limiting 
positions both of the great circle planes chosen for the growing segment 
[e, La — e] and the decreasing one [— e, La + e] as e — > 0. (We assume that 
the definition segment for the curve A may be infinitesimally extended in 
both directions.) Formally, this choice is a plane spanned by tAo, ^Ao + ^Ai- 
For such a closure, 71 + 7o = tt and Eq. (^) becomes 

WrA = —TwA — - mod 1, 

while Eq. (|2^ does not change (as usual, by Sa+g is meant the area swept 
out by the closed curve A + G defined as above). If, in addition, UAo + nAi = 
0, then a plane spanned by t^o and hAo has to be chosen. 

Remark 2. In the above consideration we have used the ribbon gener- 
ally based on the principal normal (cf. |^]), though any other continuous 
ribbon may be taken to obtain essentially the same formula for writhe (of 
course, the specific expressions for the angles 71 and 70 should be appropri- 
ately modified). 

In particular, the fractional part of writhe equals the twist of a special 
ribbon such that its generating normal vectors at the ends have they prop- 
erty that they could be transformed to each other by a parallel transport 
along the closing geodesies (i.e. 71 + 70 = mod 2it). 

5.2 A broken curve 

The above approach may be naturally extended to a sequence of disjoint 

n 

segments. Let A = ^ A^'^^ be a set of n continuous fragments. Each A^'^^ is 

i=l 
(i) (i) 

oriented so that Aq and A\ be its initial and terminal points, respectively. 
Based on A^^\ a ribbon R^^ may be built as it was done for a single piece of 
curve. We also construct n additional pieces that tie the ending point of the 
i-th fragment to the initial point of the consequent one. We identify formally 
the point A^^^^^ with A^^^ to make the entire curve closed. The connecting 
parts are built in exactly the same way as the closure of a single segment 
was made in the previous subsection. Thus, we can repeat our arguments 
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to obtain 

n 1 " 

WrA = -J2'^w,--Y,il?+ll^^) mod 1, (31) 



27r 

i=l i=l 

where Twi is the twist of the ribbon R^^ , the angles ''y[^^ and 7g'^ are deter- 
mined by their trigonometric functions 

cos^W = _^ ^i sinjP = ^^4L2^=, (32) 

^AO ^Al J Y V ^0 ^1 

cos^W = ^ ^^o-'^i , sin7W = — =E ^o"^^i (33) 



and t^^^,n^)^, 

'^'■^ Frenet frames at the beginning {v — 0) and at the 
end {v = 1) of the i-th segment (z = 1, 2, . . . , n, j = 1 + {i mod n)). 

Note that the value of writhe generally depends on both the order of 
fragments and the orientations along them. 

5.3 A non-smoothly closed loop 

A particular case takes place when the segment A forms a non-smoothly 
closed shape. This means discontinuity of the tangent vector at the initial 
point and the tangent indicatrix is not closed. The whole procedure de- 
scribed above may be well applied to such a loop though one complication 
appears: the resulting closed curve to which the basic Calugareanu- White- 
Fuller formula is to be applied has a self-intersection point at the beginning 
of the loop considered. Generally, the writhe is not defined for such shapes. 



However, on the one hand, it was shown in Section 2.3 that the GauB in- 
tegral exists unless the tangent at the loop starting point directs exactly 
opposite to the end tangent. On the other hand, under the same limitation, 
we can restrict ourselves to consideration of two limiting curves approaching 
the self-intersection shape from two different sides. As it is well known, the 
writhe jumps by 2 as a curve crosses itself [20|. Thus, the fractional part 
of the writhe is not affected by self-intersection and may be computed by 



Eq. (26) or Eq. (28) in the same way as for the open segment. 
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6 Adding the writhe 

The aim of this section is to obtain an equation that expresses the writhe of 
a segment concatenated from two or more shorter segments as a function of 
the writhes of those smaller elements. 



6.1 A special case of a closed tangent indicatrix 

We start with a consideration of an open or closed segment A = r^(s) : 
[0, L^] — > M'^ that has a (smoothly) closed tangent indicatrix A on S^. Let 
D £ Ahe the point corresponding to the starting tangent tyi(O) and a point 
E £ A represents some other value of s = si : < si < La- We draw a 
geodesic G to tie the points E and D (oriented from E to D). Denote the 
two subsegments Ai = rA{s) : [0, si] and A2 = tca{s) : [si, L^] M?. 

Thus, Ai -\- A2 = A. The writhe of the initial segment A is, according to 
Eq. (ID, 

WrA = — - 1 mod 2. 
27r 

Apply now Eq. (|2^) to the both parts of A: 



WrA, = mod 1, WrA, = mod 1. 

(Recall that the sign "— " in —G denotes its reversed orientation.) 
Clearly, Sai+g + Sa2-G = Sa and therefore 

WrA = WrAi + Wr_A2 mod 1. (34) 

We can also reformulate this addition rule by expressing the writhe as a 



function of the twisting numbers. Then, Eqs. (24), (pq) imply 



WrA = —TwA mod 1, 
Wr^i = -TwA^ - ^^'^^^ mod 1, WrAj = -Twa2 + J" ^° mod 1. 

Here, 71 is an angle from n/i(si) to to-go, the initial normal vector of G. 
Respectively, 70 is an angle from nci, the terminal normal vector of G, to 
nA(0). 

The twist Tw is additive, hence Twa = Twa^ + Twa2 ^^'^ again come 
to Eq. IM). 
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6.2 Another special case: two fragments 

If the segment A forms a non-closed tangent indicatrix A £ S'^, then the 
addition rule is more complex. Let F £ A, F ^ D he the point corresponding 
to t^(L^). In addition to G = ED, we draw two more geodesies: H = FE 
and K = FD (Fig. ^). Applying Eq. ( p6| ) to the entire AAi and to its parts 
A2 yields 

WrA = ^^ mod 1, 
27r 

Wm, = %±^ mod 1, WrA, = ^^^ mod 1. (35) 




Figure 5: The tangent indicatrix is divided by point E into two fragments 
Ai and A2; G, H, and K are geodesies. 

Taking into account the additivity of area, we have 

Sa+k = Sai+g + Sa2+h + Sadef mod in. (36) 

The triangle DEF is formed by geodesies and its signed area may be calcu- 
lated as 

Sadef = (Jdef{^DEF + ZEFD + ZFDE - vr), (37) 



E.L.Starostin. On the writhe of non-closed curves 



23 



where aoEF = sign((tAo x tAm) • t^i), t^m = tA(si)- 

The angles of the triangle can be found knowing the vectors of its vertices 
D, E, F (t AO, tAm and tAi, respectively): 

cos K — cos G cos H 

cos /LDLr = , 

sin G sin H 

cos G — cos H cos K , , 

cos AEFD = . = , (38) 

sin H sin K 

cos H — cos K cos G 

cos AFDE = = = , 

sin K sin G 

where 



cos G = tAO • tA-m, cos H = tAm ' tAl, COS K = tAl- tAO, 

sinG = |tAo X tAmI, sin# = |tAm X tAil, sin^ = |tAi X tAol- 

Coming back to the addition rule for the writhe, we make use of Eq. (p^. 
The resulting equation then takes the form 

Wta = WrAi + WrA2 + mod 1. (39) 

By means of Eqs. (|37|), (|3^ , the area term in the last equation may be 
computed explicitly. 

Note that ii E E K, then the area term in Eq. (^) is zero. 



6.3 The writhe of a growing curve 

Though Eq. (^) is valid modulo 1, it can be used for computation of the 
exact value of the writhe of a curve which can be considered as being incre- 
mentally elongated. Suppose we add to the curve Ai a fragment A2 which 
is short enough such that jWrAj H — ^^r^\ ^ 1- Assume also that the exact 
value of WrA^ is known. In the absence of self-intersections, the writhe is 
a continuous function of the arclength and we can apply Eq. (^) omitting 
mod 1 to get the exact value of the writhe of the extended curve. Starting 
with a short fragment and reiterating the above procedure, it is possible to 
compute the writhe of an arbitrarily long curve in n steps where n is the 
number of increments. 
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6.4 The general case 

We now come to the general case. Let ^ = r^(s),0 < s < be a non- 
closed segment with the tangent indicatrix A £ S^. Let {si},i = 1, . . . ,N, 
be the values of the arclength such that < si < S2 < ■ ■ ■ < Si < Sj+i < 
. . . < Sat < La- For every Si we can compute the tangent tAi = tyi(si) and 
denote t^o = ^^(0) and t^_jv+i = ^a{La)- Let £ A, i = 0, . . . , N + 1, 
correspond to tAi- 

Choose some j, 1 < j < iV, and apply the addition rule to the two 
segments EqEj and EjEn+V- 

Wro,^+i = Wroj + Wr,- ^+1 + mod 1, 

where 5oj-,Ar+i is the signed area of the geodesic triangle AEqEjEn^i. If 
J > 1, we can further apply the same rule to compute W^qj as a sum of 
the writhes of subsegments. The same is also possible to do with Wrj^jv+i, 
if j < N. After iterating this procedure as far as possible (until all the 
Wrj^j+i, i = 0,...,N, are present in the summation formula) we finally 
come to 

Wro,^+i = V Wn,+i + mod 1, (40) 

^ — ' ZTT 

i=0 

where 5*0,1, ...,Ar,Ar+i is the signed area of the spherical polygon formed by 
geodesies Eqi, £'12, • • -,Ejj^i, . . .,En,n+i, En+i,o (by Eij is denoted the 
geodesic EiEj, i ^ j)- This area may be computed as a sum 

N 

So^i^,,,^N,N+i = '^S{AEQEiEi+i) mod An. 
1=1 

By using Eq. ( |37| ) the area of each spherical triangle can be presented in the 
form 

S{AEQEiEi^i) = cTo,j,j+i(Z£'o-E'j-E'i+i + ^EiEi+iE^ + ZEi^iEoEi — tt), 

o"o,i,j+i = sign((tAo X tAi) ■ t^i+i), 

and each angle can be expressed as a function of tAi (see Eq. (^)). 
Another way to represent the area of a spherical + 2-gon is 

N+l 

So,i,...,N,N+i = ^ cFi-i^i^i+ii'^Ei^iEiEi+i - vr) mod 2-k, 

i=0 

identifying E^i = E^+i and En +2 = Eq- 
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7 The writhe of a polygonal line 

In this section we consider a curve consisting of a number of straight line 
segments. Such a curve can serve as a model of a self-avoiding walk [10| or 
as a skeleton description of a linear polymer Q. 

Though the polygonal line has a discontinuous tangent, it may be con- 
sidered as a limit of a smooth curve which differs from the polygonal line by 
arbitrarily small planar arcs in the vicinities of every point of discontinuity. 
The consequitive tangent directions are then connected by geodesies on 



[16|. It can also describe the path of a scattered light beam [pO| ]. 

Conversely, an arbitrary smooth curve may be approximated by a polyg- 
onal line. In one can find an estimate of the difference between the 
writhing numbers of a closed smooth curve and a polygonal curve inscribed 
within. 

Several methods have been proposed for computation of the writhe of 



closed polygons |32, 33]. Here we give an exact expression for the fractional 
part of the writhe for an arbitrary polygonal line. The formula does not 
involve the double summation. 

Let A = {Pi}, z = 0, . . . , + 1, be a sequence of points in M'^. For every 

interval {Pm, Pm+i), a tangent vector is = PmPm+i/\PmPm+i\, m = 

> > 

0, . . . ,N. We assume that PmPm+i H PnPn+i = v, if t„ x = 0, Vm ^ 
n, m,n = 0, . . . , N. The tangent indicatrix ^ G is a line consisting of 
fragments of geodesies EmEm+i, m = 0, . . . , N — 1, Em = tm r\ S'^ £ A. 
According to the above developed approach, we add the closing geodesic 
EnEq] it corresponds to an additional (A^ + l)-th fragment of A such that 

^N+l = to- 

In order to make use of the formula for writhe, we have to calculate the 
signed area of a spherical polygon V = {Eq . . . E]\f}. It may be represented 
as a union of triangles 

N-l 

V= [j AEoEiE,+i. 

i=l 

Each triangle AE^EiEi^i has a signed spherical area 



'0,i,i+l 



crQ,i,i+iS{AEoEiEi+i) mod in, 



where S{AEoEiEi+i) = ao,i,i+i + ai,i+i,o + ai+i,o,i — is a non-negative area 
of a spherical triangle with the angles akim and cTkim = sign((tfc x t/) • 



E.L.Starostin. On the writhe of non-closed curves 



26 



The signed area of the polygon V is 



N-l 

Sv = 2. 'S'o,i,i+i 



mod 47r. 



i=l 



The angles akim are functions of tangents such that 




We now are able to rewrite Eq. ( p6| ) in this particular case: 



y^TA = ^ X] '^o,j,j+i mod 1. 



(41) 



i=l 



This expression is valid for any curve in that generates a tangent 
indicatrix containing only geodesic fragments. Also, it seems to be useful 
as an approximation while calculating the writhe of a smooth curve which 
can be properly discretized. The last means that the discretization should 
be better performed on its tangent indicatrix rather than on the curve in 

itself. Conventional methods of integration may be applied to compute 
an approximate value of the area swept out on S^. 

8 The writhe and the Gaufi integral 

The writhe of the smooth closed curve may be expressed as the double 
integral (Eq. @)). It is evident that the writhe for an open segment of the 
length L as defined above can no more be computed as the Gaufi integral 
over this segment though, in most cases, the double integral itself is also 
well defined for smooth non-closed curves. 

Our aim here is to obtain a formula connecting the both values: on 
the one hand, the writhe that relates to the difference between the linking 
number and the twisting for the ribbon built with the geodesic closure and, 
on the other hand, simply the double integral taken over the open segment. 

8.1 An open curve and its closure 

Consider an open smooth non-self-intersecting curve A = r(s) : [0, L] — > R'^. 
We assume here that the tangent vectors t(s) = r'(s) are neither parallel nor 
antiparallel at the ends: t(0) / ±t(L) (we will examine these cases later). 



E.L.Starostin. On the writhe of non-closed curves 



27 



We extend the curve A with two straight hne segments: B : rsisi) = r(L) + 
sit(L), si G [0,/] and C : rc(s2) = r(0) + S2t(0), S2 G [-/,0]. Note that the 
both segments have the same length I. Now connect the end points of B 
and C with the straight hne segment Di : r£)(^) = (1 — S^)rB{l) + £,rc{—l) = 
(1 - + ^t(L)) + ^(r(0) - lt{0)). The direction of Di is determined by 

its tangent 



toil) 



\rci-l)-rB{l)\- 



Now let the lengths of B and C increase to infinity and compute the 
limiting orientation of the tangent t £> : 

Thus, we see that, in the limit I ^ oo, Dqo lies in the plane defined by 
the initial and end tangents of the segment A. 

In the case when t(0) = ibt(L), we can also attach two straight line 
segments. If t(0) = — t(L), then all straight lines connecting these segments 
belong to the same plane defined by t(0) and r(L) — r(0). The case t(0) = 
t(L) requires a special consideration. 

What we have now is a closed circuit A + B + Di + C. It is smooth except 
for two points at the beginning and at the end of Di. We modify B and Di 
in the small vicinity of where they join themselves together by introducing 
a planar curvilinear segment E with the tangent varying from t(L) to t/j. 
All three segments involved belong to the same plane spanned by t(L) and 
ti). We can assume that the length of E does not depend on /. The length 
of the shortened segment B^, is decreased to be 

The similar procedure may be carried out to smoothen the join of the 
segments Di and C. The new planar curvilinear segment F belongs to 
the plane spanned by t(0) and tn- The length of F is the same for every 
I; without loss of generality, we assume that the length of the shortened 
segment equals /*, as well. 

We have come to the smooth closed curve A + + + + F + C* . 
We are interested in the limiting case when — > oo. The tangent indicatrix 
of the initial curve A is closed then by a geodesic corresponding to the 
limiting curve E^o + -Doo + -^oo • This follows from the construction of these 
curves and from Eq. (42). The limiting curve -Bqo + -^'oo + -Doo + -Foo may be 
considered as an implementation of the first part G of the closure constructed 



in Section |5Jj. Thus, the writhe of the open segment A may be computed 



E.L.Starostin. On the writhe of non-closed curves 



28 



as the writhe of the Hmiting closed curve A + B^o + E^o + D^o + -Foo + Co 



and its fractional part satisfies Eq. ( ^q) and Eq. (28). 

However, for the smooth closed curve A+Boo+Eoc + Doo+Foo + Coc, the 
writhe may be obtained by the double integral formula independently. Since 
the circuit consists of 6 parts, we are to consider all the pairs of them as 
they are involved in the double integration. For the brevity, we will denote 
an integral over a pair of curves P and Q by (P, Q). Clearly, (P, Q) is the 
same as {Q,P). 

Before proceeding with this, we obtain some simple estimate of the value 
of the double integral 

C2 Ci 





where 



j j Iwri^y 1,(^2) daida2, 



(ri((Ti) - r2((T2)) • (ti(cri) X 12(0-2)) 



|ri(o-i) - r2(o-2)r 

The integral I2 is taken over two smooth curves ri{ai),ai £ [0,£i] and 
r2(cr2),cr2 G [0,£2]- Let A = min |ri((Ji) - r2(cT2)| > 0. Then 

<. ...... <. IIj,;^^^ s («) 





Eq. ( ^31 ) implies that lim I2 = for any two curves of finite length. If 

A—* CO 

one of the curves has its length of order A or less, i.e., Ci = 0{A),i = 1, 2, 
and the other is of the finite length C^-i, then the integral I2 vanishes as 
A 00, too. 

We now come back to the integral over the pairs of curves. The inte- 
grals over {B,,B,), {B,,E), {B,,D,), {E,E), {E,D,), {D,,D,), (Z)*,F), 
(D*, C*), (F, F), {F, Gk), (C*, C*) equal zero because the integrand vanishes 
for coplanar curves. If t(0) = —t{L), then the whole closure is planar and the 
integrals {B^:,F), {E,F), {E,C^:) also vanish for every l^,. If t(0) / ±t(L), 
then the length of is of order for large Therefore, by applying 
Eq. (^3|), we obtain that the integrals {B^,F), {E,F), {E,C^) as well as 
{A,E), (A,D^:), {A,F) all approach zero as 00. We denote the remain- 
ing integrals that can be non-zero as follows: 

yVy = ^ j j Iwris,s) dsds, 

A A 
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Swi = ^ j j Iwris,si) dsids, SW2 = -^ j j Iwris,S2) ds2ds, 

Boo ^ Coo 

Sq = ^ j j />Vr(si,S2) dsids2. 

We call them "the wrying", "the swirl" and "the squint", respectively. 
Thus, the writhe of A may be represented as a sum 

y\;r = Wy + Swi + Sw2 + Sq. (44) 

The first summand Wy is simply the double integral over the open seg- 
ment under consideration. Therefore, Eq. ( ^4| ) provides a connection be- 
tween this integral and the writhe. 

Let us now examine the case of parallel tangents t(0) = t(L). The two 
attached segments B and C have opposite directions then and, instead of 
the straight line -D, we connect them by a circular arc D joining B and C at 
the same distance I from the ends of A. The arc D lies in the plane spanned 
by the vectors r(L) — r(0) and t(0) and its length is of order / for large /. 
Again, the smoothening curves E and F can be constructed in the similar 
way as it is done in the regular case. Thus, we obtain the smooth planar 
closure of the curve A. 

After letting the lengths of B^, and C=k go to infinity and analysing the 
double integral components in the expression for the writhe we come to the 
same Eq. (p^) with the right-hand terms defined as above. 

It may occur that the ray B ox C intersects the curve A. Then, generally, 
the writhe of the whole closed curve A+B+E+D+F +C is not determined. 
The situation is the same as for a non-smoothly closed loop (Section [5.31 ). 
In the generic case, when the tangents in the point of the intersection are 
neither coincident nor of opposite direction, the fractional part of the writhe 
still can be found by the examination of the two limiting positions of the 
curves in the vicinity of the intersection point. Since the writhe jumps by 2 
as the curve goes through itself, the half-sum of the writhes for those curves 
may be taken as the value of the writhe. The same approach may be applied 
to another singular case when the rays B and C cross each other. Moreover, 
the constraint of non-self-intersection of the initial open fragment A my also 
be weakened in the similar fashion. 

Next we are going to clarify the structure of the integrals Sw and Sq. 
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8.2 The swirl 

Consider 

L oo 

27: j j |r(s)-r(L)-sit(L)|3 



L 



27r 

Here we denoted 

oo 

ds\ 



)- j {Yi{s)-{t{s)xti)isn,{s) ds. (45) 



ISw{s) 



[(R-(R-ti)ti)2 + (R.ti-si)2]2 



and R = R(s) = r(s) - r(L), ti = t(L). 
We can represent Isw{s) as 



CXD 

/ds 
^, 6 = 6(s) = R-ti, = a2(s) = (R-6ti)2, 



[a2 + (6-si)2 
and carry out the integration to get 

si — h 



Isw(s) 



a^y/a^ + {b-siy 



+ 62(7^2 + 52 _ 5) • 



(46) 



Substituting a(s),5(s) into Eq. (|4q ) and the result further into Eq. 
yields 

1 R(s) • (t(s) X ti) 



'^^^"2^70 |R(.)|(|R(.)|-R(.)-ti)^'- ^^^^ 

Let us introduce the spherical coordinate system with the origin at the 
point r(L) and let the z-axis be directed along the ray B (Fig. |^. Then 
R(s) = (pcos^/^cosc/), pcos^/^sin0, /JsinV'), ti = (0,0, 1), and p = p{s), (p = 
<j){s), ip = ip{s) are the functions describing the curve A. 

In these coordinates, Eq. (|47|) takes the form 



Swi = — I 4>'{1 + sin V')fis. (48) 
2vr J 
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Figure 6: The spherical coordinates p,(p,ip. The z-axis is directed along the 
tangent t(L) at the end point of the curve r(s). 
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Note that the swirl is zero if (j){s) = const, i.e., if the whole curve A 
is planar. The swirl also vanishes when ip{s) = — ^ which means that the 
curve A is a straight line continuation of the ray B. 

It is natural that the swirl is scale- invariant: it does not depend explicitly 
on how remote are the points of the curve from the ray. 

The second integral Sw2 over the ray Coo has the same structure. 



8.3 The squint 

Now consider the integral Sq over the two rays i?oo and Coo ■ It is convenient 
to introduce the special Cartesian coordinates with the origin at the point 
r(0) and the x-axis directed along r(L) — r(0) (Fig. 0). Let y-axis lie in the 
plane of the ray B^o and the z-axis be chosen such that the whole coordinate 
system is right-handed. Denote hy € [0, vr], the angle from the x-axis to 
the direction of t(L). The orientation of the ray Coo is defined by two angles: 
ip,ij £ [0, vr], between the x-axis and t(0) and 9,9 £ [0, 27r], between the xy- 
plane and the plane spanned by the x-axis and t(0). In this coordinate 
system we may represent the both rays as follows: 

rB(si) = (5 + Sicos0,sisin(/),O), g = |r(L) - r(0)|, 
ts = (cos 0, sin 0, 0), si £ [0,oo], 
rc(s2) = (s2 cos ^/^, S2 sin V' cos 0, S2 sin V' sin 0), 
tc = (cosV',sin'(/'cos0,sinV'sin6'), S2 G [— oo,0]. 

We wish to compute the integral 

00 

^ 1 f f {rB{si)-rc{s2))-{tB{si)xtc{s2)) ^ , 

Sq= — / -— 13 dsids2 = 

2vr y J \rB{si) - rc{s2)\-^ 

-00 

00 

1 ^. \ f f rnisi) - rc{s2) , , 

= *^ X / / ] — j-^ dsids2 = 

27r J J |rB(si) - rc(s2)r 

-00 

00 



= — sin -0 sin (/) sin / / ds-idso, 

2vr J (s2 + 2psi + g2)i 

where p = p{s2) = gcos<j) — S2{cos ip cos cj) + sinV'sin(/)cos0), = q^{s2) 
s| — 2gs2 cos ij) + g^ . 
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Figure 7: The angles (p,9,ip determine the orientation of the tangents t(0) 
and t(L) at the ends of the curve segment. 
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It is easy to perform the first integration: 







oq = — sm yj sm sm u 






= — sin w sm (psmB / — ^ 

27r ^ J q{p + q) 

— oo 



— sin Tp sin sin 9 



27r J ^Js2 — 2gs2 cos ip + g'^ 



-oo 



dS2 



5 cos (j) — S2 (cos COS + sin -0 sin (j) cos 0) + ~ ^552 cos tp + 

The last integral can also be done and the result may be presented as 
an algebraic formula which docs not depend on g, naturally. However, the 
derivation of the final expression involves a complicated algebra and, instead, 
we prefer to obtain Sq in a different way. 

We can consider the both rays and the straight line connecting them as an 
(infinite) polygonal line with three links. The writhe of this line, as defined 
in the previous section, is exactly equal to Sq. Then the squint is essentially 
proportional to the signed area of the spherical triangle constituted by the 
geodesies that tie the vertices corresponding to the vectors t(0), r(L) — r(0), 
and t(L). 

The triangle has its two sides equal to and cf) and the angle between 
them 0. By the side cosine theorem for spherical triangles, we find the third 
side X from 

cos X = cos tp cos (?!) + sin ^ sin cf) cos 6 

and the signed area of the triangle can be calculated by the generalized 
Heron's formula 



S T,-(f> T.-ip S-x 

5 = 4i^ arctan \ / tan — tan — - — tan tan — - — , 

V 2 2 2 2 ' 

where u = sign((r(0) — r(L)) • (t(0) x t(L))) = sign(sin'0sin0sin^) and 

E = i(0 + V^ + x). 

Then the squint is Sq = 
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9 Example: helical shapes 

In this section we shall be dealing with a regular circular helix: r(s) = 
(cos as, sin as, - a^s), < a < I, < s < L, s is the arc coordinate 
and L the length of the segment. The limiting values of the parameter a 
correspond to a straight line (a = 0) and a circle of the unit radius (a = 1). 
The curve is periodic with the period T = The tangent, normal and 
binormal vectors are 

t = (—a sin as, ocos as, \/ 1 — a^), n = (— cos as, — sin as, 0), 

b = (a/ 1 — a^ sin as, —\J\ — cos as, a), 

and K = a^, r = a^J\ — a^ are the curvature and the torsion, respectively. 
The tangent indicatrix of the helix is a circular arc on S'^ of radius a. 



9.1 The writhe of an arbitrary segment of the helix 

We now illustrate the above developed approach by applying Eq. ([2^ ) to the 
helix. We need to calculate the twist first. The normal n(s) is a well-defined 
continuous vector function and 

L 

Tw{L) = Twf{L) = ^ [ Tds = ^a^/l-a^. (49) 
zvr J zvr 



The second step is to calculate the angles 71 and 79. To this purpose, 
we find the tangent, normal and binormal vectors at the beginning (s = 0) 
and at the end (s = L) points of the segment: 

to = (0,a,A/l-a2), no = (-1,0,0), bo = (0, -x/T^, a). 



ti = (—a sin aL, a cos aL, yl — a?), ni = (— cos aL, — sin aL, 0), 



bi = (a/ 1 — a^ sinaL, — \/l — a^ cos aL, a). 
It is possible now to make use of Eqs. (p9|),(|30|) to obtain 



cr cos ^ 

cos 71 = cos 70 = , (50) 





a2 sin2 ^ 




- a2 sin ^ 



V J- Uj oi±i n 

sin 71 = sin 70 = , (51) 



1 - a2 sin2 ^ 
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where a = signsin^. Eqs. (|50|),(|5l|) imply 71 = 70 and tan 71 = tan 70 = 
-VI -a2 tan^^. 

After substitution of these angles and the expression for Tw from Eq. ( |49| ) 
into Eq. (pSD we obtain the fractional part of the writhe 



Wr{L) 



1 

TT 



arctan ( \/l 



a? sin 



aL aL\ 
— , cos — 
2 ' 2 I 



L 



mod 1. 



We denote by z = arctan(2;,y) a function such that sin 2; 
-vr < z < TT. 



■■ , cos z 



Since the writhe of the helix should be continuous, we may rewrite the 
last equation in the form 



>Vr(L) 



1 



vr 



arctan 



Vl-a2 



aL 



sm — , cos — I av \ — a? 

2 ' 2 ) 2 ^ 



+ 



+2roundf^j, (52) 



where round(x) is a function that gives an integer nearest to x. 

The writhe Wr as a function of the arclength L normalized on the period 
T is presented in Fig. ^ for three different values of the parameter a. 

We may use Eq. (^) to compute the writhe for the particular lengths 
L = ^ [m is the number of half-periods): 



Wri — 



\ m , J — 
)=y(l-vT 



(53) 



(cf. |2g,g]). 

In the limiting case a ^ 0, Eq. (52) produces Wr — > as it should be for 
a straight line. Another limiting value of the parameter is 1: as a — > 1, the 
helix approaches a circle r(s) = (cos s, sins, 0), covered L/27r times. Then 
Eq. ( |52[ ) gives Wr round (^)- 



9.2 The writhe by the double integral 

We wish to make use of the double integral formula (see Eq. (^) in order to 
compute the writhe of the infinite helix by an alternative means. We take a 
helix A with 2M turns, M > 1 being integer, and close it as described above 
with the planar curve C. (We could equally take an odd number of turns, but 
the even number facilitates the algebraic manipulations.) We can build C as 
a union of three fragments: C = Ci + C2 + C3 . All the fragments lie in the 
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L/T 

Figure 8: The writhe Wr of the circular hehx as a function of the arclength 
L normahzed on the period T = In/a, for three different values of the 
parameter a. 
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plane x = 1. Ci starts at the end of the hehx in the point (1, 0, MTV 1 — a^) 
and it has the initial tangent {0,a,^/Y^^a^). Denote its length by Li and 
let it end at the point {l,g,MT\/l — a^), g = const, g < Li, such that the 
whole curve is convex and it lies in the half-plane x = 1, z > MT^Jl — a?. 
The end tangent vector is (0, 0, —1). 

The fragment C2 is simply a straight line r2(cr) : (1,^,-0"), a G [—MT 
VT^^, MTVl - a2] . 

The last part C3 is similar to Ci. It begins in the point (I, g, —MT 
\/l — a^) with the tangent (0, 0, —1) and ends exactly in the beginning of the 
helix (1, 0, — MT\/l — a^) having the tangent (0, a, Vl — a^)- The convex 
curve C3 has the length L3 and it belongs to the half-plane x = 1, z < 
-MT^l-a"^. 

It makes no sense to try to get the writhe for the infinitely lengthy helix. 
Instead, we consider the writhe per period defined as 

^^(^) = J"fL^-i / / I^ris,a)dsda, (54) 

A+CA+C 

where 

, . ^ (r(g)-r(a))-(t(g) xt(a)) 

Iwr{s,a) = -— --|3 . (55) 

\r[s) — r[a)[^ 

In Eq. (^), we imply that the curves Ci and C3 remain unchanged 
as M varies and only the length of the straight line fragment C2 increases 
appropriately with M growing. Then it may be shown that the right-hand 
side of Eq. (|5^) depends on the shape of neither Ci nor C3. 

The double integral in Eq. ( |54|) may be represented as a sum 



A+C A+C A A A 

(we took into account that j j = J j\ i = 1,2,3, and / / = 0, Vi,j = 

A Ci Ci A Ci Cj 

1, 2, 3, since lyvr = at every point). 

Now consider the mixed integral over A and Ci. We split it into two 
parts 

(Af-l)TLi MT Li 

I\\;rdads= J J Iw/rdads+ J J I\\;rdads. (56) 

A Ci -MT (M-l)T 
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The second integral in Eq. ( pq ) is taken over two finite fragments of a smooth 
curve and its value is finite, too, and it does not depend on M. As to the 
first integral in Eq. (^6[), we can estimate its numerator 

\{r{s) - r^{a)) • m X t,{a))\ < 



< V (Li + 2)2 + (Li + \/l - a^{MT - s))2 < 
< 2(Li + 1) + - a^{MT - s) 

and the denominator 

|r(s) - ri{a))f > (1 - a^)^MT - sf. 

Then 



(M-l)T Li 



Jyvr dads 



-MT 



(Af-l)T 



-MT 



Li 



2(Li + 1) + Vr^^(Mr - s) 
(l-a2)f(Mr-s)3 

2AfT , 

2(Li + 1) + yr^^^ 



(l-a2) 



(l-a2 









2Mr y 



< 



< 



(l-a2)f 



2^2 



Thus, we have shown that the mixed integral over A and Ci is bounded 
uniformly with respect to M. The same property for the integral over A 
and C3 can be proved similarly. We now deal with the third mixed integral 
over A and C2. It may be explicitly written as 



MT MTVT^ 



1 — COS as — g sin as 



(2 + 3^2 _ 2 cos as — 2g sin as + — a'^s + (t)2 



dads. 



-mt_mtVT^ 

Since we are interested in obtaining the limit of the integral as M ^ co. 
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we can consider 



MT 



l2 = —a J (1 — cos as — g sin as) X 

-MT 



da 



(2 + g'^ — 2 cos as — 2g sin as + — a^s + o")^) 2 
It is easy to check that 



ds. 



di 



^, P>0. 



Therefore, 



MT 



-2a 



1 — cos as — g sin as 
2 + g"^ — 2 cos as — 2g sin as 



ds. 



-MT 



The last integral may be explicitly computed 

1 — cos as — g sin as 



2 -\- g"^ — 2 cos as — 2g sin as 



ds 



s 1 

2 ~ a 



as 



arctan ( — ((4 + y ) tan — — 2g) 
9^ ^ 



vr , / as 
— round — 
a \2'K 



It is periodic in s with the period T and, since 



s vr 1 / as \ 

round — 

2 a \2tx) 



< 



IT 

2a 



and 



as 



arctan — ((4 + (7 ) tan 2g) 



5" 



< 



TT 



2a' 



we get I/2I < 47r. 

Summing up, we have proved that only double integration over the infi- 
nite helix A may be retained in Eq. (54) 



Wr{T) = lim 



1 1 



MT MT 



M^oo 2MT 4tt 



In>r{s, a)dsda. 



(57) 



-MT-MT 
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For the problem at hand, the integrand /yvr may be reduced to 



/wr(s,cr) = a\/l — a? x 
2sin(f(s -0-)) [2sin(f (s - a)) - a(s - a) cos(f (s - a))] 



[4sin2(f(s-(T)) + (l-a2)(s-a 



^ 3 
,2 2 



It depends only on the difference of the arc coordinates due to the translation 
invariance along the helix curve. It is natural to make the change of variables 
(s, cr) — > (u, w), u = s — a, v = s + a and to consider 

; , , ? . X / r 2sin(fn) r2sin(fn) - aucosffu)] 

[4sin2(fn) + (l-a2)u2]2 

For any fragment of the helix which length equals the period, the writhe 
is the same and it may be found as 

MT MT oo 
nV(r) = Jim^ ^ • i- j dv j iyvr{u)du = ^ j iy^r{u)du. (58) 

-MT -MT -oo 



The justification of why Eq. (57) is equivalent to Eq. (pq ) follows from 
Lemma 2 (see Appendix B) applied for (/? = 0. 
Consider an integral 



^ / T / \ 7 T / \ sva.w[sva.w — w cosw 
h= Jo{w)dw, Jo{w) = 



(sin^ w + Xw'^) 



where A = const, A > is a parameter; we put A = 
For small w, 



3 ' 
2 



1 \w\ a 



Mw) = ' '3 + 0{w') = -\w\ + 0{w'). (59) 
3(1 + A)2 3 

The integral Ii may be represented as a sum 

— e e 00 00 e 

Ii = j Jo{w)dw + J Jo{w)dw + J Jo{w)dw = 2 J Jo{w)dw + J jQ{w)dw, 

—00 — e € e — e 

where e > is small. 
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In view of Eq. (59) 



lim / Jo{w)dw = 



and we get 



OO 

2 lim / Jn{w)dw = 2 lim 

e^O J e^O 



W 



a/ sin^ w + Xw"^ 



Returning to Eq. (|58D and making the identification w = |n, we obtain 



Wr{T) 



a(l - Vl - a?) Wr(T) 



27r 



T 



(60) 



The last expression agrees with the value of the writhe for m = 2 (Eq. (|53D). 

We remark that the computation of the writhe (and twist) of the helix 
is a favourite example of an application of the first Fuller theorem (e.g., see 
|0|, H, m) though the author has not come across the direct usage of the 
double integral formula in the literature. It is evident, that the last approach 
is much more tedious, but it seems to be instructive to see it at work. To 
complete the picture, we show in the next section how to calculate the same 
quantity by means of the Fuller second theorem. 



9.3 The writhe by the FuUer second theorem 

We consider two curves. The first one is a segment of a straight line 
Aq : tq = (1, 0, Vl — a^s) and the second is one period of the helix Ai : 
ri = (cos as, sin as, Vl — a^s), the same as above. The common parameter 
s varies from to T = — (note that it is not the arclength coordinate for 
Ao). 

The Fuller second theorem is only applicable to closed curves, so we have 
to close both our segments. We shall do it in accordance with the procedure 
described above for periodic curves. Namely, we close ^-nd Ai each with 
a curve lying in the plane x = 1 and consisting of three fragments. The 
closure of Aq contains the parts: 

1) Coi, a semicircle of radius i?o > with the tangent toi = (0, — sin(/), 
cos(/>), < (p < TT,(j) is a new parameter; 
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2) Co2, a segment of a straight hne parallel to Aq: ro2 = (1, 2Rq, — Vl — 
{T-e)),9e[0,T]; 

3) Cos, a semicircle of radius Rq with the tangent toa = (0, — sin cos 
vr < (/> < 27r. 

The closing curve for the helix is built as follows: 

1) Cii, a semicircle of radius i?i > vrVl — with the tangent tn = 
(0,— sin((/) +(/)o), cos(</) + 0o)), < (/> < vr, (/>o is defined by — sin</)o = 
a, cos</>o = \/l — a^; 

2) Ci2, a segment of a straight line: ri2 = (l,a(a^ — 1)^ + 2RiVl — a^, 
-(1 - a2)l0 + ^^T^^ - 2Ria), 9 e [0, T]; 

3) Ci3, a semicircle of radius R2 = Ri — 7r\/l — with the tangent ti3 = 
(0, -sin(0 + 0o),cos((/) + (/>o)),vr < < 27r. 

At the joining points the tangent vector is not differentiable, but we may 
always make the closing curves smoother by small variations in the vicinities 
of the joining points. These modifications influence the result in no way. 

The integrand in Eq. (^) vanishes on the straight line segments of the 
closure and the rest two integrals over the semicircles cancel each other 
because toi(</>) = — to3(<^ + vr) and tn^cp) = —ti^^cp + tt) for < < vr. 

We now may forget about the closing segments and take the only integral 
over the initial range of parameter s E [0, T]. The straightforward calcu- 
lations yield the expression for the integrand function ^^^-^ ^2 = a(l — 

Vl — a^) and we get again Eq. (60) for the writhe per period. In contrast 
to the previous section where we dealt with the infinite helix, we now have 
obtained the same formula by considering only one helical period. 



9.4 A double helix 

The procedure of the calculation of the writhe of a helix is very similar to 
that for a double helical shape. The helix is assumed to be regular, circular 
and closed at both ends in the same manner as it is described above for 
open curves. The difference is that each of the two closing curves joins two 
different helices. 

In this section we show how to obtain the limiting value of the writhe per 
turn when the integer number of turns tends to infinity. We include into our 
consideration the case when the double helix is not symmetric with respect 
to the central axis (as it takes place in the B-form of DNA, for example). 
Each of the two helical curves will be called a strand. 

The first strand is described as the helix Dq : ri(si) = (cososi, sinasi, 
— a^si), < a < 1, So < si < Si, si is the arc coordinate and 
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L = Si — So the length of the segment (cf. the beginning of Section ^). 
The second strand may be obtained from Dq by rotation through a (con- 
stant) offset angle around the central axis: —D2 : r2(s2) = (cos(as2 + 
2(p),sm{as2 + 2tp), Vl — 0^*2), Sq < S2 < Si, S2 is the arclength. (Note 
that it is convenient to choose the same orientation for the parametrization 
of both curves so that —D2 is oriented against -Do-) 

Let the double helix have 2M turns, M integer (i.e. L = 2MT, T = 
and Si = MT = —Sq. It is closed with two loops Di and D3 of finite lengths 
Li and L3, respectively. We assume that the loop Di is entirely contained 
in t he cyli nder xj + yf < 2, zi > \/l - a^MT and D3 G {xj + yi < 2, 2:3 < 
— Vl — a'^MT}. The whole closed double helix may be represented as the 
sum D = Do + Di + D2 + D3. 

We wish to compute the writhe per period 

"^^<^' = Jsi. 2LTra -ill 

D D 

where I-\/\>r{s,a) is as in Eq. ( |55|). 

The double integral in Eq. ( |6l|) may be represented as the sum 



D D ^'^-^bi Dj 

Neither the integral / J nor J J depends on M, hence they do not affect 

Di Di D3 D3 

the limit in Eq. (|6l|) . Consider the mixed integral over Dq and Di. We can 
split it similarly as it was done for the single helix (cf. Eq. (p^) 

(M-l)TLi MT Li 

I\\;r dads = J J Iy\^rdads+ j j ly^rdcrds. (62) 

Do Di -MT (M-l)T 



The last double integral in Eq. ( |6^ ) exists and does not depend on M. As 
to the first integral in the right-hand side of Eq. (|6^), we can estimate its 
numerator 

|(ro(s)-ri(a))-(to(s) xti(a))| < 



< Y (1 + 2)2 + (Li/2 + Vl - a2(MT - s))2 < 3 + Li/2 + y^l - a2(Mr - s) 
and the denominator 

|ro(s)-ri(a))|3>(l-a2)l(Mr-s)3. 
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Then it is easy to check that 

(M-l)T Li 

Iy\!r dads < 



-MT 



Li + 6 



(l-a2 



4r2 



+ 



We can conclude that the mixed integral over Dq and Di is bounded 
uniformly with respect to M. The same is true for the pairs {D-^^Dq), 
{Di,D2), {02,0^) and, of course, for (1)1,1)3). Now we see that the only 
pairs that count in Eq. ( |6l| ) are those involving both strands. Note that the 
integral / / = / / was already estimated as / J (see Section |9.2|) . 
Do Do D2 D2 A A 

The last integral to be computed is / / = — / / . We apply Lemma 2 

Do D2 Do -D2 

(see Appendix B) to obtain its limiting value which, according to Eq. 
will be 



^ / Iwr,ip{u) du 



w 



4vr ^/^jp^w^^^pyT^^ 



27r 



a 
'2^' 



au l-a^ 
w = -, A = — ^. 

Note that the result does not depend on ip. 

Adding the above value to Eq. (^ ) (which is actually the limiting value 
of the double integral over the same helix), we finally obtain the writhe per 
period for the double helix 



Wr{T) 



i\/r 



27r 



(63) 



It is necessary to clarify that the above formula gives the writhe normalized 
on the whole length of both strands. Note also the negative sign of the 
writhe for the right-handed (a > 0) double helix. 

The symmetric case of 2(/? = vr was considered in |34], where essentially 
the same formula for writhe was found. 

We remark that the writhe for n-strand helical shapes (n > 3) may be 
easily computed on the basis of the results derived for the single and double 
helices. 
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9.5 The writhe of the double hehx of arbitrary length 

Here we continue dealing with the non-symmetric double helix, but now we 
consider the general case when the helix does not need to have an integer 
number of turns. The length of each strand is the same and denoted by L. 

We apply the equation for writhe that was found earlier for the sequence 
of disjoint fragments (Section 5.2). In the double helix case, we have only 
two pieces: Ai = Dq and A2 = —D2, they are defined in the previous 
section, but now we set Sq = and Si = L. The twist of each helix is the 
same 

L 



Twi = Tw2 = ——a^/l — a 



2n ^ 



(cf. Eq. (H)). 

Now compute the vectors of the Frenet bases at all four ending points 

aW aW .(2) ^(2). 

t^') = (0,0,^1-02), n^') = (-1,0,0), b5J = (0,-\/l-a2,a), 
(1) 



t^] = (— asinaL, acosaL, \/l — a2), n^] = (— cos oL, — sin aL, 0), 

b^] = (\/r^^^sinaL, — vT^^^cosaL, a). 



t^Q = (asin(aL + 299), — acos(aL + 2(^), — v 1 — a?), 
n^Q = (— cos(aL + 2(p), — sin(aL + 2(/?), 0), 
b^Q = (— v'T^-^sin(aL + 2ip), \J \ — cos(aL + 299), —a), 
*Ai ~ (a sin 299, —a cos 2(/?, — — a2), n^] = (— cos 2^9, — sin 2(/?, 0), 

(— Vl — a2 sin2(/?, y \ — cos 2^9, —a). 



-(2) 



It is easy to see that the angles 7q*\ 7^*^ i = 1,2, computed with the 
help of Eqs. (p^ ), (|33|), satisfy the following equations: 



which imply 



(1) (1) (2) (2) 

cos 7i = — cos 7q = cos 7i = — cos 7q , 

(1) . (1) . (2) . (2) 

sm ^\ = sm 7q = sm 7^ = sm 7q , 



7^ ^ + lo^ + 1? + if^ = mod 27r 



and according to Eq. (^) we have 



L 



Wr = a\/ 1 — a2 mod 1. 
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This expression is in compliance with the writhe per length for an integer 
number of turns (Eq. (|63|)). If we think about the double helix as a contin- 
uously growing structure, then it is evident that the last equation gives not 
only the fractional part of Wr, but its exact value as function of one strand 
length. 

Let h be the length of the axis of the double helix, h = L\/l — a?. Then 

Wr = 

TT 

We see that the growing double helix delivers an example of a family of 
curves A{h), parametrized with the continuous parameter /i, such that the 
writhe is a linear function of the length. Clearly, the writhe per (double 
helix) length is constant. Also, note that the writhe does not depend on the 
offset angle (p which controls the mutual location of the strands (in terms of 
DNA we may reformulate the last observation as an invariance property of 
writhe with respect to the widths of the minor or major grooves). 



10 Conclusion 

We have been concerned with the generalization of the notion of the writhe 
for an arbitrary space curve and with obtaining effective formulas for its 
computation. We have analyzed various conditions on the position and 
orientation of the ends of the fragment in space starting with a periodic 
curve and finishing at the most common case. 

In all cases, an explicit construction of a closed ribbon was carried out. 
This makes it possible an application of the Calugareanu- White-Fuller equa- 
tion to the determination of the writhe. A rule for the calculation of the 
writhe for a segment which is a union of smaller parts, is obtained. We 
examined a polygonal line special case. 

A relation was established between the writhe and the Gaufi integral 
taken over the open fragment. The difference between these two quantities 
may be represented as three single integrals. 

The application of the formulas presented was demonstrated on the ex- 
amplar curves, including regular single and double helices. In particular, the 
writhe as a continuous function of arclength is defined for a regular helix. 
It was shown that a double helix of finite length, with "geodesic" closures 
at the ends, provides an example of a one-parameter family of curves that 
realizes the linear dependence of the writhe on the length and the writhe is 
invariant with respect to the value of the offset between the strands. 
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In Appendix A, a new derivation of the formula connecting the writhe 
and the twist for a non-closed ribbon with one continuous edge is given. 

11 Appendix A: A formula for a non-closed ribbon 

The formula, we are about to deduce, relates to a ribbon based on a smooth 
closed curve but generated by a vector function which has a discontinuity. 
Such a ribbon may serve as a model to a nicked circular DNA, one strand 
of which is cleaved. 

We will consider a closed smooth non-self-intersecting space curve A = 
r(s) : [0, L] of class C^. We assume for the sake of simplicity that 

the parameter s is the arclength. The smooth closure implies r(0) = r(L), 
to = r'(0) = r'(L) = ti, r"(0) = r"{L). Let the curve A be equipped with 
the continuous vector function u(s) : [0, L] such that u(s) • r'(s) = 

0,Vs G [0,L]. 




Figure 9: A non-closed ribbon based on the smooth closed curve A. 

In the general case, u(0) ^ u(L) and we can define an angle a between 
these two vectors, measured up to modulo 27r. This angle is a characteristic 
of discontinuity by which the ribbon (r, u) fails to close (Fig. ^). 
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Our aim here is to obtain an equation connecting the writhe of the curve 
A, its twist and the angle a. 

Consider the tangent indicatrix A of the curve A on S^. It is a closed 
curve. Choose a great circle plane containing to = ti. It has at least one 
common point with the tangent indicatrix. We look for a plane having at 
least two more intersections with A. Such a plane almost always exists. 
Indeed, take a great circle which is tangent to A at the initial point with 
the tangent tg. Generically, there is an alternative: either the great circle 
intersects A (the first case) or not (the second one). In the first case, it is 
enough to rotate the great circle plane about to through a small angle to 
provide two additional points of intersections that bifurcate from the initial 
one. Consider the second case. The whole tangent indicatrix of a closed 
curve may not lie in only one of the two semispheres (see Theorem 1 in 
[p5| ). Then the closedness of A implies at least two crossings in addition to 
the common initial point. 

We should distinguish the singular case when the tangent indicatrix is 
itself a great circle. It corresponds to the planar space curve. We shall 
return to this case later. 

For a regular configuration, we can choose two additional points P2 ■ 
r{s2) and P3 : r(s3) for which to+pt2 + qts = 0, p,q are real and t2 ^ r'('S2)j 
t3 = r'(s3) (it is easy to show that the points may be always taken such 
that t2 + t3 / 0). 

The starting point and P2 and P3 divide the curve A into 3 parts, namely, 
A = A1 + A2 + A3, Ai : r(s), s £ [0, sa], A2 : r{s),s £ [s2, S3], A3 : r{s),s £ 

[S3,L]. 

We now are ready to construct another curve by inserting 3 new segments 
of straight lines between the parts. The new curve B consists of 6 fragments 

Di : r{0) + ato, a£ [0,ei]; 

Bi : r(fT - ei) + eito, a £ [ei, ei + S2]; 

D2 : r(s2) + eito + (fi - ei - S2)t2, cr G [ei + S2, ei + S2 + £2]; 

B2 : r{a - ei - €2) + eito + e2t2, a £ [ei + €2 + S2, ei + €2 + S3]; 

D3 ■■ r(s3) + eito + e2t2 + ((T-ei-e2-S3)t3, a £ [ei + £2 + 53, ei + e2 + S3 + e3]; 

S3 : r((T-ei -£2-63) + eito + e2t2 + e3t3, a£ [ei + e2 + e3 + S3, ei + e2 + e3 + L]. 

The compound curve B{ei, €2, €3) = Di + Bi + D2 + B2 + D3 + B3 is closed 
if eito + e2t2 + €313 = 0. We put €2 = pei and €3 = qei and choose e > 
such that B{ei) would have no self-intersections for Vei, ei £ [0, e]. Thus, we 
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have constructed an one-parameter family of closed pulled-out curves B{ei), 
ei G [0,e], B{0)=A. 

The next step is to define ribbons (i?(ei),U), Vei G (0, e]. For the frag- 
ments Bi, i = 1,2, 3, we keep the same vectors U = u as for the respective 
parts Ai, i = 1,2, 3. For the fragments D2 and D3, U = u(s2) = const and 
U = u(s3) = const, respectively. For the remained fragment Di, we define 
U(o-), a G [0,ei] such that U(0) = u(L) and U(ei) = u(0): 

a a 
U((TJ = u(L) cos — a + (to X u{L)) sm — a. 

To find the angle a between u(0) and u(L) we have the relationship 

u(0) = u(L) cos a + (to X u(L)) sin a. 

The ribbon (i?(ei),U) is continuous and the Calugareanu- White- Fuller 
formula can be applied to get 

£^B(ei) = WrB(ei) + TwB{e{). (64) 

Since we have chosen e small enough to exclude self-intersections of the 
curve (more accurately, we have to require an absence of self-crossings for 
the whole ribbon), the right-hand side of Eq. ( [6^ ) does not depend on ei. 

By the construction, the tangent indicatrix of B{e\) is the same for 
any ei (namely, it is A because -B(O) = A). Therefore, due to the Fuller 
first theorem (and taking into account the absence of self-intersections), 
WrB(ei) = Wr^. 

The twist of B[e\) may be represented as a sum 

TwB{ei) = TwDi(ei) + Twbi + TwD2{e.i) + Twb2 + Twdz{<^i) + Twbs- 

Clearly, Tw£)2{ei) = and Twosiei) = 0, because the vector U is constant 
on these straight-line fragments. Further, Twbi + Twb2 + '^ss = '^w^, the 
twist of the initial non-closed ribbon (which is well defined). 

The twist of the remaining part Di may be readily computed 

a 

and it does not depend on ei. 

We now can rewrite Eq. (|6^ ) as 

CkB = WrA + TwA + -^. 

zvr 
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We may interpret the last equation as a definition of the linking number for 



the non-closed ribbon A (cf. the last equation in |2S] and Eq. (7.2) in pO|| ): 



a 

CkA = WrA + TwA, CkA = Cks - —■ 

ZTT 

In conclusion, we explain how one has to treat the degenerated case of a 
planar curve A. Clearly, two insertions are sufficient then. Indeed, we can 
choose the second point on A as diametrically opposed to the starting one, 
i.e., to + t2 = 0. In other words, we can proceed the same way as in the 
regular case after putting p = \ and g = 0. It is evident, that the result will 
not change. 

12 Appendix B: Estimation of integrals 

In this section two lemmas are proved that are useful for computation of the 
Gaufi integral for regular helices. 
Lemma 1. Let 

sin((w — 9)[s\n{w — 6) — w cos{w — 9)] 
" [sin\w - 6) + Xw^f/^ ' 

where A, 9 are parameters and A > 0. 
Then the integral 

oo 

-^6»oo = j wJe{w)dw 



exists. 

Consider an integral 

wo 

Ie{wo) = j wJQ{w)dw 



and integrate it by parts using the equality 

top 

Jg{w)dw 







yjsm^{w - 9) + Xw'^ 
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Then 



leiwo) = = - iiwo), 



sin^(t(;o - 9) + 



i{wo) = J 



wo 

w 



^ \/ sin^(u) — 0) + Xw"^ 



dw. 



Now we estimate /(luo) for wq > 0. On the one hand, 

dw wq 



Wo 

i{wo) < / ^ = ^ 







and on the other, 



wo 

i{wo) > J - 





Therefore, we have 



ivdtv Y^l + A»'q — 1 



"° "° < /,(-o) < 



sm'^{wo-e) + Xwl ^ Y^sin2(w;o - 0) + Xw'^ 

yi + Aw^ - 1 
A 

Letting wq go to the infinity in the last inequality leads to 

< leoo = lim Ie{wo) < Y 

which proves the statement of the lemma. 
Lemma 2. Let 

2sin(|(s — a) — ip) [2sin(|(s — a) — ^p) — a{s — cr) cos(|(s — a) — (p)] 



[4sin2(|(s - a) - (^) + (1 - a'^){s - af] ^ 

Then 

MT MT oo 
-MT -MT -oo 

(65) 
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where 



I 2sm(^n — 99) r2sin(f n — 99) — aucos(f u — 

^y^r,^\u) =a\/\- f -. 

[4sin2(§u-^) + (l-a2)u2]2 

First of all, note that /wr,(/p(S) f) = iv^r-Lp^—s^—a) (and Iwr,ip{u) 
Iy\;r-^{—u), as well). Next, we make the change of variables (s,cr) 
{u,v), u = s — a, V = s + a to obtain 



MT MT 



-MT -MT 



MT MT 



Iwr,<fi{u)dudv— (66) 



-MT -MT 



MT u 





2MT 2MT-U 



{iwr,^{u) + iwr-<fiiu))dvdu 



MT 



Due to Lemma 1, 

MT u 



MT 



1 



lim , 

M^oo MT 




1 



Iwr,±<p{u)dvdu = lim — — / uIyyr,±Ju)du = 





It is also evident that 
lim 



2MT 2MT- 



1 



M-^oo MT 



Iwr,±ip{u)dvdu 



MT 
2MT 



1 



lim 

M^oo MT 



{2MT - u)Iv^r,±^{u)du = 0, 



MT 

-2\ 



because Iwr,±ip{u) = 0{u~ ) for u oo. 

Therefore, we may consider only the limit of the first pair of the integrals 
in the right-hand side of Eq. ( |66| ) which will be 



MT MT 



1 1 



lim 

M^oo 2MT 4tt 



Iwr,ip{u)dudv = ^ I Iwr,ip{u)du. 



-MT -MT 
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